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The effects of multi-impurity on the entanglement of anisotropic Heisenberg ring XXZ under a 
homogeneous magnetic field have been studied. The impurities make the equal pairwise entangle- 
ment in a ring compete with each other so that the pairwise entanglement exhibits oscillation. If 
the impurities are of larger couplings, both the critical temperature and pairwise entanglement can 
be improved. 
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INTRODUCTION 



Entanglement is not only the fabulous feature of quan- 
tum mechanics but also very important to the quantum 
information processing (QIP).[^1 In the studies of quan- 
tum entanglement, solid state system with Heisenberg 
model interaction is the simple and applicable candidates 
for the realization of quantum information. Therefore, 
there are many works focusing mainly on the different 
kinds of Heisenberg models [^^^^1 such as spin ring etc.. 

The impurities often exist in solid system and plays 
a very obvious and important part in condensed mat- 
ter physics. As a candidate of QIP, solid system with 
impurity is also one of our important study object. 
In the previous researches, the impurity effects on the 
quantum entanglement have been studied in a three- 
spin system [22,23] g^,^^ ^ large spin systems under zero 
temperature. [^'^1 However, in these works, they have just 
studied single impurity. 

In this paper, we will focus on studying the effects of 
multi- impurity on the pairwise thermal entanglement in 
a ring chain. We find the impurities make the equal pair- 
wise entanglement in a ring compete with each other. If 
impurities are of large couplings, the critical temperature 
and the pairwise entanglement which coupled to the im- 
purities can be improved. Our studying results not only 
provide a standard to judge impurities but also provide a 
way to enhance entanglement and critical temperature. 




FIG. 1: Two configurations of spin ring when the impurities 
are non-nearest neighbors, (a): qubit ring formed with 10 
qubits. The 4th and 6th are two identical impurities, (b): 
qubit ring formed with 10 qubits. The 4th, 6th and 8th are 
three identical impurities. 
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II. NON-NEAREST NEIGHBORING IMPURITY 
EFFECT 



(1) 



Firstly, we investigate the multi-impurity effect when 
the impurities are non-nearest neighbors shown as the 
Fig. 1. In the two figures, square represents impurity 
qubit and round stands for normal qubit. 

For the case of Fig. la, the Hamiltonian can be written 

as 



where (erf, a^, erf) are the vector of Pauli matrices; J 
and Jz are the real coupling coefhcients of arbitrary near- 
est neighboring two qubits. We restrict the B > along 
z direction and N + 1 = 1. We choose the parameters B, 
J, Jz and T are dimensionless and assume the coupling 
coefficients between normal qubit and impurity one has 
the relation 
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(a) 



J = a* J,J, = a* J, 



(2) 



where a characterizes the relative strength of the extra 
couphng between the impurity and its nearest neighbor- 
ing qubitsJ^^l For the case of Fig. lb, one can write it 
easily following Eq.(l). We do not give it here any more. 

As we know, for a system in equilibrium at tempera- 
ture T, the density operator is p = [l/Z] cxpi—H/ksT), 
where Z — rr[exp(— /f/fcsT)] is the partition function 
and ks is Boltzman's constant. For simplicity, we write 
ks — 1- The value of entanglement between two qubits 
can be measured by Concurrence C which is written as 



C = max(0, 2 max — A^) 



(3) 
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[25,26,27,28] which Ai is the square roots of the eigenval- 
ues of the matrix 



(4) 



where p is the density matrix and the symbol * stands for 
complex conjugate. The Concurrence can be calculated 
no matter whether p is pure or mixed. In the following, 
we just take the pairwise entanglement into account. We 
will trace over the qubits and study the reduced density 
matrix of the two qubits which we are interested in. 
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FIG. 2: Nearest neighboring concurrences as a function of 
a for the two-impurity model (4th and 6th are two identical 
impurities) (a) and the three-impurity model (4th, 6th and 
8th are three identical impurities) (b). T=l, B=0.4, J=l, 
J^=0.65. 

Now, we review the difference between an ideal ring 
chain and an open chain. For an ideal ring chain, every 
qubit is of the same position with the others so that any 
pairwise qubits are of the same amount of entanglement. 
But for an ideal open chain, pairwise entanglement is re- 
lated to the position of the qubits and exhibit oscillations 
due to the breaking of the symmetries. ^^^^ 



FIG. 3: Nearest neighboring concurrences as a function of 
i for the three-impurity model (4th, 6th and 8th are three 
identical impurities) q=0.1 (a), q=2 (b). T=l, B=0.4, J=l, 
J^=0.65. 



Because here we study multi-impurity, we can not ob- 
tain analytic expression of the system. We will directly 
numerical calculate and plot entanglement. In Fig. 2, we 
plot the pairwise entanglement as a function of a, corre- 
sponding to Fig. la and Fig. lb, respectively. For both of 
the two cases, in the regions of far away from impurities, 
the entanglement, for example C12 and Cioi, are slightly 
affected by the various values of a. Within the impu- 
rities regions, we observe that there is the almost same 
threshold value of a, after which a qubit and its nearest 
impurity start entanghng such as C34,C45 etc.. In Fig. 
3, we show clearly that the pairwise entanglement versus 
site i. If a is small shown in Fig. 3a, the case equal to 
cutting at 4th and 8th, thus the chain 9-10-1-2-3 is sim- 
ilar to the open chain^^^l while the part 4-5-6-7-8 chain 
have no entanglement because of the weak couplings. If 
a > 1 such as a = 2 shown in Fig. 3b, we still can cut the 
chain into two parts because J' > J. Within the pure 
regions, entanglement will compete while in containing 
impurity part pairwise entanglement still compete each 
other. 

Fig. 4 shows the influence of temperature and the val- 
ues of a on the entanglement in three-impurity model. 
From this figure, we can judge again that the second and 
the third qubit are pure qubits while the third and the 
forth contain one impurity. Usually, it is difficult to ad- 
just the coupling coefficients, which means we will meet 
with difficult if we directly use the behavior of Fig. 3 
to judge which one is impurity. But we still can do it by 
measu ring the Concurrence changing with temperature 
(Refs |29l [3(| proposed that Concurrence can be mea- 
sured), because changing the temperature is very easy. 
On the other hand, we find that a can effectively en- 
hance the Concurrence and critical temperature if a > 1 
which is show in Fig. 3 and 4 clearly. By introduc- 
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FIG. 4: Nearest neighboring concurrence C23 and C34 versus 
a and T for the three-impurity model (4th, 6th and 8th are 
three identical impurities). B = 0.4, J = 1,J^ = 0.65. 




j=4,6 
i=l 

with 

j" =P*J,j'^ =p*J,, (6) 

where /3 characterizes the relative strength of the ex- 
tra coupling between the two nearest neighbor impuri- 
ties and J , still has the relation of Eq.(2). Similarly, 
one can write the Hamiltonian corresponding to Fig. 5b. 
We plot the pairwise Concurrence near the two-nearest- 
impurity qubits area as a function of /? which is shown 
in Fig. 6. From this figure, we can see easily that the 
nearest neighboring impurities coupling only affect the 
nearest two-impurity and the others which couple with 
the impurities. For example, in Fig. 6a, the nearest neigh- 
bor impurities C56 has a threshold value of /?, affected 
by /? heavily while C45 also will decrease as a results of 
the competition between neighbor qubits. For the case 
of Fig. 5b, although we have more impurities, the nearest 



FIG. 5; Two configurations of spin ring with nearest-neighbor 
impurity, (a): qubit ring formed with 10 qubits. The 5th and 
6th arc two identical impurities, (b): qubit ring formed with 
10 qubits. The 4th, 7th and 8th are three identical impurities. 



ing impurities with large coupling, one can also improve 
critical temperature and entanglement. Therefore, our 
studying not only provide a stand to judge impurity but 
also exhibit a way to enhance entanglement and critical 
temperature. 



III. NEAREST NEIGHBORING IMPURITY 
EFFECT 
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In this section, wc study the nearest neighboring im- 
purity effect on entanglement. We study the rings with 
a structure of Fig. 5. According to the Fig. 5a, 

the Hamiltonian is 



FIG. 6: Nearest neighboring concurrences versus P for the 
two-nearest-impurity model (5th and 6th arc two identical 
impurities) (a) and the three-nearest-impurity model (4th, 
7th and 8th are three identical impurities) (b). B=0.4, J=l, 
J^=0.65, a=0.8. 
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neighbor coupling only affect entanglement of themselves 
C78 and that coupling with the impurities Cgr, Cgg; and 
all the others pairwise entanglement almost can not be 
affected. 



IV. CONCLUSION 

In conclusion, for a Heisenberg XXZ ring under a ho- 
mogeneous magnetic field, we have studied entanglement 
in two- impurity and three-impurity under the two case of 
non-nearest-impurity and nearest-impurity. We find that 
the introducing of impurities make the originally equal 
pairwise entanglement compete with each other. For the 
weak and strong a, we can cut the ring chain into dif- 
ferent open chain and then use the open chain property 



to explain the competition. For the case with nearest 
neighbor qubits , the change of the relative coupling /? 
can only affect the qubits which couple to the impurities. 
If introducing impurity with large a and /?, the pairwise 
entanglement, which couple with the impurities directly, 
can be enhanced and critical temperature also will be 
improved. 

This work was supported by Natural Science Founda- 
tion of China under Grant No. 10575017 and Natural 
Science Foundation of Liaoning Province of China under 
No. 20031073. 



V. REFERENCES 



Bennett C H and DiVincenze D P Nature 2000 404 247 [17 
Nielsen M A Phys. Rev. A 2001 63 022114 [18^ 
Kamta G L and Starace A F Phys. Rev. Lett. 2002 88 
107901 [19 
O'Connor K M and Wootters W K Phys. Rev. A 2001 
63 052302 [20 
Dong L L and Man S L Chin. Phys. 2002 11 383 [21 
Yeo Y Phys. Rev. A 2002 66 062312 

Sun J R, Wei Y N and Pu F K Chin. Phys. 1995 4 542 [22 
Zhang G F, Li S S and Liang J Q Opt. Commun. 2005 [23 

245 457 

Zhou L, Song H S, Guo Y Q and Li C Phys. Rev. A 2003 [24 

68 024301 [25 
Wang X G Phys. Rev. A 2001 64 012313 

Gunlycke D, Kendon V M and Vedral V 2001 Phys. Rev. [26 
A 64 042302 [27 
Arnoson M C, Bose S and Vedral V Phys. Rev. Lett. 2001 [28 

87 017901 

Canosa N and Rossignoli R Phys. Rev. A 2004 69 052306 [29 
Khveshchenko D V Phys. Rev. B 2003 68 193307 
Zhang Y M and Xu B W Chin. Phys. 1995 4 842 [30 
Zhang T, Hui X Q and Yue R H Acta Phys. Sin. 2004 53 
2755 



Shang Y M and Yao K L Chin. Phys. 1998 7 864 

Sun Y, Chen Y G and Chen H Phys. Rev. A 2003 68 

044301 

Shao Y Z, Lan T and Lin G M Acta Phys. Sin. 2001 50 

948 

Dong Z H and Feng S P Chin. Phys. 1998 7 348 

Hui X Q, Chen W X, Liu Q and Yue R H Acta Phys. 

Sin. 2006 55 3026 

Fu H, Solomon A I and Wang X J. Phys. A 2002 35 4293 
Xi X Q, Chen W X, Hao S R and Yue R H Phys. Lett. 
A 2002 297 291 

Wang X G Phys. Rev. E 2004 69 066118 

Bennett C H, DiVinccnzo D P, Smolin ,] A and Wootters 

W K Pliys. Rev. A 1996 54 3824 

Wootters W K Phys. Rev. Lett. 1998 80 2254 

Hill S and Wootters W K Phys. Rev. Lett. 1997 78 5022 

Antcnoodo C and Souza A M C J. of Opt. B:Quantum 

Scmiclassical Opt. 2003 5 73 

Sancho J M G and Huelga S F Phys. Rev. A 2000 61 
042303 

Davidovich L "Entanglement As An Observable" 2006 
reported in Texas A&M University 



